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ABSTRACT 


H« iMresent a tachnlque for the meaauraaent of magnetic heliclty from 
values of the two point magnetic field correlation matrix under the assumption 
of spatial homogeneity. Knowledge of a single scalar function of apace, 
derivable from the correlation matrix, suffices to determine the magnetic 
helicity. We illustrate the technique by reporting the first measurement of 
the magnetic helicity of the solar wind. 


A oonsiderabl* body of theoretical plasma physios literature over the 
last twenty.five years has emphasized the importance of magnetic helioity, 
defined as 


z J’A*jBd*x 


( 1 ) 


where JB and ^ are the magnetic field and the vector potential, respectively. 
This integral extends ov«r all field containing regions, and jk is subject to 
the gauge condition s 0. may also be defined "per unit voluse" so that 
\ z The magnetic helicity measures the departure of a turbulent 

magnetic field from mirror symaetry, or equivalently, the degree of topologi- 
cal linlcage of magnetic flux tubes'. 

Vbltjer* noticed Uiat under fairly gmeral assumptions H is an integral 
invariant of the incompressible one fluid ideal magnetohydrodynamio (MHD) 
equations. He then developed a variational formulation for calculating MHD 
equilibria* which Tbylor* and Montgomery ^ al.* applied to MHD dynamics in a 
conducting cylinder. Their "relaxation theory" has been used in the MHD 
theory of Reversed Field Pinch plasms confinement devices*. 

Turbulence theory has also addressed the dynamical role of magnetic 
helicity. Frisch^ al.*, in analog!> to two-dimensional turbulence theory*'*, 
conjectured that an inverse cascade of magnetic helioity may be characteristic 
of three-dimensional MHD flows. In steady state, with helioity and energy 

supplied at a constant rate to intermediate wavenumbers, an MHD system would 
then consist of two inertial ranges: an inverse cascade of magnetic helicity 

to large scales, and a direct cascade of energy to amall scales. A class of 
selective decay hypotheses has bem discussed* *>• which states that the large 
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Reynolds number, non steady dynaaios, of a two rugged invariant system are 
oharaoterized by the rapid decay of the direct transfer quantity and (reserva- 
tion of the baok transfer quantity. 

The Qonjeotured implications for three-dimensional NHD are that: 1), the 
ratio of energy to magnetio helioity decreases in time; 2), in the limit of 
ex^emely high Rsynolds mabers the ratio approaches a limit defined by the 
geometry; and 3), the net magnetic helicity .resides primarily in the largest 
scales of the system, regardless of its scale in the initial conditions. 

Our understanding of the dynamical importance of and the range of 
validity of the relaxation, inverse cascade and selective decay theories is 
limited by the lack of interplay between theory and experiiient. Nb are not 
aware of a single direct measurement of magnetic helicity.or its speotrua. 
Thus, it is desirable to develop procedures for obtaining magnetio helioity 
from experiments, m this letter we resent a straightforward {w^ocedure for 
evaluating magnetio helioity from values of the two point magnetic field 
correlation matrix under the assumption that the statistical properties 

are spatially homogeneous. 

Vli begin with the definition 

Rij(0 a (2) 

where the brackets denote an appropriate average over a statistical ensemble. 
The assumption of "wesk" homogeneity renders a function only of the 

spatial separations^; thus ®ij satisfies the solenoi- 

dal condition 


( 3 ) 



Mh«re suBBation is implied by repeated indices. As discussed in Ref. 11, 
Rij(^) may always be additively decomposed into a symmetric proper tensor 
Tij(^) with even spatial parity and an antisymmetric pseudotensor 
odd q)atial parity. 

The transformation properties of a homogeneous, solenoidal correlation 
matrix R^^ uniquely determine the form of the pseudotensor can be 

seen by use of the method of isotropic tensors introduced by Robertson and 
others**. By using all available vectors and the isotropic forms and 
'iJR, exhaustively list the linearly independent tensor forms 

available for inclusion in R^^. The required manipulatiois are easier to 
perform in Fburier apace employing Uie energy speotrtss tensor 

SijCk) . (^)- /dr .-’i'l («) 

If no rotational or refleotional symmetries arn assumed, the solenoidal 

constraint leads to 31 k-spaoe tensor forms at the onset**. These depeiwl on Jk 

and two independent principal axis unit vectors. The Fourier space pseudo- 

tensor, odd under k must be antisysmetric in its indices; similarly, the 

prcper tensor is even and symmetric. Simple algebraic manipulations ^w that 

only three pseudctensors can be linearly independent. The solenoidal 

constraint eliminates two of these leaving one form, c, , k G(k), where G is a 

1 jm m •* 

scalar function of satisfying G(Jc) 9 . Additional sysnetries imposed 

on th«« system can only modify the way in which G depends on but cuinot 
introduce additional functions into the form of the k-spaoe pseudotensor. 
Likewise in configuration space, is derivable from a single scalar 

function «(r), and has the form 
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• ‘ija Jr 

B 

1h« faot that dapanda only on tha gradiant of a single scalar function of 
r provides the basis for determining the aagnetio helioity. 

Let the ij^^ oonponant of ttM Fburiir transform of the symmetric part of 
the correlation be denoted Then, using Eq. (4), H^jCJO 


Because is the trace of <B evaluated at « 0 (Eq. 

have 


s /dk H..(k). 
m jj <i«^ 

From Eq. (6) and (7) • we nov see that 
helioity, • Furthermore , because 

syametrio part of G^c) * 

have the desired expression for vix .. 


1), we iiaediately 
(7) 

is the spectrum of magnetic 
I^(J^ depends solely on the anti- 
Finally, from Eq. (1) and (5), we 


s a*B> a 2a(^ a 0) 


( 8 ) 


The function «(r) , idiioh determines both the total magnetic helioity and 
its spectrum, may be evaluated by performing a line integral over separation 
values: 


r 

0t0 

2e(^) s 2 ^ 7#»d^ 


L ‘*^1 *lJm^jm^A^ 


(9) 
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Eq. (9) is valid provided that tha oorrelationa vani^ rapidly aa |£| ♦ 
is than given by Eq. (8)* 

1b dataraine on* Bust Fburiar tranafora Eq. (9), Mhioh requires 

exhaustive knowledge of • tor nil However, a reduced helioity speotna is 
available fV'aa knotdedge of R^^(r^,0,0). the oorrelation tensor for a aequenoe 
of ooSinear separations in the direction. The reduced energy spectrua 
tensor^* is defined as 

Sjj(k^) i ^ /dr^ e’^1*‘l R^j(r^,0,0) (10) 

s /dkgdkj Sj^jCk^.kg.kj) 


so that 

Hj[(k^) * 57 ♦(r^.0.0) e"^“l‘'l (11) 

s 2 In S230c^)/k^ 
and 

* / dk^ Hjj(k^) > 2e(r^«0,0,0) (12) 

Application of these results to fusion plasmas nay be United by strong 
inhoaogeneities present in a bouaded laboratory device. However, in the 
laboratory one can interpret the ensenble average as an average over identical, 
ly prepared "shots" of a containnent device, which nay pemit useful values of 
helioity to be extrscted. Magnetic helicity can also be neasured in space 

plasmas. Vb have used nagnetMeter data fron the Voyager 2 spacecraft in the 
solar wiixl near 2.8 astroncnicai units (AU). The results reported here are 
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froB • singlt spat of OH tows on days 95 - 97 of 1978. This is a singla 
point owasuranant it', a highly supar-illfvtaio flow, which justifias usa of tha 
MHD analogua of tha frozai flow appoximatiai to obtain s|[j(k). 

Nb avaluatad via two indapandant naans: tha BLaotaian-Tukay (BT) 

aaat laggad produot taohniqua with 24 dagraas of fraadoa, and tha Fast Fburiar 
Ttranafom (FFT) taohniqua snoothad to hava aquivalant statistical validity. 
Tha two taohniquas gava assantially idantioal rasults. lasts wars aada to 
alininata tha possibility that aliasing or laakaga affaotad tha analysis, lha 
rasults shown in Fig. 1 and 2 wara obtainad vVith tha FFT approach. 

Fig. 1 displays boUi tha raduead n^natic «aargy dansity E(k) ^ 
and |kl^(k)| (cf. Eq. 12) in units of r* (1 r » Gauss). (Ito oonvarting 
froB fraquanoy to wavaniabar, wa usad tha solar wind ^ad during thib period 
vrtiioh was 450 kn/s. in tha following discussion, tha superscript on tha 

reduced helioity is omitted. 

As has often bean reported, E(k) closely an>roxiBatas a power law (tha 
slope in this case is -1.7 i 0.1). Tto notable feature in Fig. 1 is that tha 
envelope of |kf^(k)| closely ^aoas tha sana power law. Note that there is no 
tendency for |ki^(k) |/E(k) to baooBa small at large k. m fact, kl^(k) 
oacillatas between ■ ±0.4 of its Baxiaal values throughout tha spaotriss. 
Neverthalaas, tha net helioity dansity (■ 18 G* oa) is entirely due to tha 
helioity in Uia largest scale fluctuations. Tha helioity containing length 
(2vH|^/£k|^(k)) is about 9.9 * 10 cm compared with tha similarly defined 
energy containing length (2vE/ZkE(k) ) , which is ■ 5 » 10** ca. The correla- 
tion length** (/R^^dr/E) for this data sat is 2.5 x 10** cm. 

Tha high degree of correlation and antioorralation between E(k) «id I^(k) 
is illustrated in Fig. 2 where k®^^(k) and k®^^H (k) are plotted against both 

B 

frequency and wavammibar. Because of tha linear scale, Fig. ? emphasizes tha 
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hif^t fr«qu«noy behavior of f^Oc). Thaaa raaults auigaat that ^ aolar Mind 
during thia period is rioh in he]i.ioal struoturea; exhibiting a aignifiorat 
twist to the magnetic field at very large aoalea. Large oluipa of poaitive 
and Mgative tMlioity are also fotaid at all smaller aoalea, but with an 
average value oloae to zero. A detailed analyais of the aolmr wind magnetio 
field struoture at several looationa in the helioaphere will be presented in a 
more complete paper. 

Ml thank the Voyager magnetometer team for help «h 1 cooperation in the 
data analysis, H. H. Mish for aid with the nvaierioal analyais, and D. 
Montgomery for many atimulating conversationa. Ihe participation cf C. Siith 
waa supported in part by National Aeronautics and Space Administration grant 
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PZGURE CAPTIONS 


Flf. 1. Th« r«duc«d aagnttlo tn«rgy d«n«lty E(k) md th« r«duo«d htlioity 

apcctruB |kH^(k)| (In «n«rgy units) of th« solar wind at 2.8 AU. Tha solar 

-12 

wind valoeity as 450 ka/a, and tha total fluotuatlon anargy is 4.6 >* 10 
arga/oa*. E(k) has a powar law alopa of k"'*^ * For clarity, not all 

valuaa of I^(k) *>** plottad at high fraquanoias. 

Fig. 2. E(k) md H^(k) plottad on a linaar aoala for tha data Aown in Fig. 
1. Tha top traca is k^^^E(k); tha bottoa traca is k®^^H^(k). 
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